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Abstract 

Flavour symmetries of Proggatt-Nielsen type can naturally reconcile the large 
quark and charged lepton mass hierarchies and the small quark mixing angles 
with the observed small neutrino mass hierarchies and their large mixing angles. 
We point out that such a flavour structure, together with the measured neutrino 
mass squared differences and mixing angles, strongly constrains yet undetermined 
parameters of the neutrino sector. Treating unknown 0(1) parameters as random 
variables, we obtain surprisingly accurate predictions for the smallest mixing 
angle, sin^ 2^13 = O.OT^qH, the smallest neutrino mass, mi = 2.2^\''[ x 10~^ eV, 
and one Majorana phase, a2i/7r = 1.0lg 2. 



1 Introduction 



It remains a theoretical challenge to explain the observed pattern of quark and lep- 
ton masses and mixings, in particular the striking differences between the quark sector 
and the neutrino sector. Promising elements of a theory of flavour are grand uni- 
fication (GUT) based on the groups SU(5), SO(IO) or Eg, supersymmetry, the see- 
saw mechanism and additional flavour symmetries [I]. A successful example is the 
Froggatt-Nielsen mechanism [2] based on spontaneously broken Abelian symmetries, 
which parametrizes quark and lepton mass ratios and mixings by powers of a small 
'hierarchy parameter' 77. The resulting structure of mass matrices also arises in com- 
pactiflcations of higher- dimensional fleld and string theories, where the parameter 1] 
is related to the location of matter flelds in the compact dimensions or to vacuum 
expectation values of moduli flelds (cf. [3J). 

In this article we consider a Froggatt-Nielsen symmetry which commutes with the 
GUT group SU(5), and which naturally explains the large z/^ — Ur mixing This 
symmetry implies a particular hierarchy pattern in the Majorana mass matrix for the 
light neutrinos. 
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which can be regarded as a key element for our analysis. The predicted Dirac and 
Majorana neutrino mass matrices are also consistent with leptogenesis [5]. Despite these 
successes, the predictive power of the Froggatt-Nielsen mechanism is rather limited due 
to unknown 0{1) coefficients in all entries of the mass matrices. For example, the 
considered model [5] can accommodate both a small as well as a large 'solar' mixing 
angle 612 [HE]. To get an idea of the range of possible predictions for a given fiavour 
structure, it is instructive to treat the 0{1) parameters as random variables [7]. 

In the following we shall employ Monte-Carlo techniques to study quantitatively the 
dependence of yet undetermined, but soon testable parameters of the neutrino sector on 
the unknown 0{1) factors of the mass matrices. Using the already measured neutrino 
masses and mixings as input, we find surprisingly sharp predictions which indicate a 
large value for the smallest mixing angle ^13 in accordance with recent results from 
T2K [8], Minos [9] and Double Chooz pH], a value for the lightest neutrino mass of 
O{10~^) eV and one Majorana phase in the mixing matrix peaked at 0:21 = vr. 
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Table 1: Froggatt-Nielsen charge assignments. From Ref. [5]. 

2 Masses and mixings in the lepton sector 

As far as orders of magnitude are concerned, the masses of quarks and charged leptons 
approximately satisfy the relations 

rrih '■ rris : ~ rrir : : rrif. 1 : rj : rf , 

with 7f ~ 1/300 for masses defined at the GUT scale. This mass hierarchy can be 
reproduced by a simple U(l) flavour symmetry. Grouping the standard model leptons 
and quarks into the SU(5) multiplets 10 = e^) and 5* = (d^, Zl), the Yukawa 

interactions take the form 

Cy = lOaO.i/, + hf^b*10,H, + h^.^5*l,H^ + ^ h^h,hS + c.c. , (3) 

where 1 = 1/^ denote the charge conjugates of right-handed neutrinos and i,j = 1...3 
are flavour indices. Note that the Yukawa matrix /i*^"^ for the right-handed neutrinos 
can always be chosen to be real and diagonal. Hu, Hd and S are the Higgs flelds 
for electroweak and B — L symmetry breaking, i.e., their vacuum expectation values 
generate the Dirac masses of quarks and leptons and the Majorana masses for the right- 
handed neutrinos, respectively. In this setup, the Yukawa couplings are determined up 
to complex 0{1) factors by assigning U(l) charges to the fermion and Higgs fields in 
Eq. ©, 

hij ~ r]^^^^^ . (4) 

With the charge assignment given in Tab. [T] the mass relations in Eq. ([2]) are re- 
produced. Additionally, perturbativity of the Yukawa couplings and constraints on 
tan /3 = {Hu) / {Ha) require < a < 1. 
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Masses 

From Eq. and Tab. [T]one obtains for the Dirac neutrino mass matrix mci and the 
Majorana mass matrix of the right-handed neutrinos M, 
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with the electroweak and B — L symmetry breaking vacuum expectation values vew 



{HuY + {Hd)"^ and vb-l = {S), respectively. In the seesaw formula 
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the dependence on the right-handed neutrino charges drops out, and one finds for the 
light neutrino mass matrix, 
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The charged lepton mass matrix is given by 
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Note that the second and third row of the matrix rrie have the same hierarchy pattern. 
This is a consequence of the same flavour charge for the second and third generation of 
leptons, which is the origin of the large neutrino mixing. Hence, diagonalizing rrie can 
a priori give a sizable contribution to the mixing in the lepton sector. 

Mixing 

The lepton mass matrices are diagonalized by bi-unitary and unitary transformations, 
respectively. 
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with V^Vi = V^Vr = WU = 1. From Vl and U one obtains the leptonic mixing matrix 
f^PMNS = ^L^' which is parametrized as [TT] 
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with Cij = cos 6ij and Sij = sin Oij. Since the hght neutrinos are Majorana fermions, all 
three phases are physical. 

In the following we study the impact of the unspecified 0{1) factors in the lep- 
ton mass matrices on the various parameters of the neutrino sector by using a Monte 
Carlo method, taking present knowledge on neutrino masses and mixings into account. 
Naively, one might expect large uncertainties in the predictions for the observables of 
the neutrino sector obtained in this setup. For instance, the neutrino mass matrix is 
calculated by multiplying three matrices, in which each entry comes with an unspeci- 
fied 0{1) factor, cf. Eq. ([6]). However, carrying out the analysis described below and 
calculating the 68% confidence intervals, we find that in many cases our results are 
sharply peaked, yielding a higher precision than only an order-of-magnitude estimate. 

3 Random variables 

Monte- Carlo study 

The unknown 0{1) coefficients of the Yukawa matrices h^'^\ h^'^^ and h^^^ are constrained 
by the experimental data on neutrino masses and mixings, with the 3a confidence ranges 
given by 

2.07 X 10"^ eV^ < lAm^^^l < 2.75 x 10"^ eV^ , 

7.05 X lO^^eV^ < AmL < 8.34 x 10"^ eV^ , 

- sol- ^^^^ 

0.75 < sin2(2^i2) < 0.93, 
0.88 < sin2(2^23) < 1 ■ 

In the following we explicitly do not use the current bound on the smallest mixing angle 
(^^13 < 0.21 at 3(7 [H]). This allows us to demonstrate that nearly all values we obtain 
for 9i3 automatically obey the experimental bound, cf. Fig. [H 

In a numerical Monte-Carlo study we generate random numbers to model the 39 real 
parameters of the three mass matrices^ The absolute values are taken to be uniformly 
distributed in [10"^/^ 10^/^] on a logarithmic scale. The phases in h^'^^ and h^'''^ are 

^Nine complex 0(1) factors in each h^'^'^ and h^''\ as well as three real C(l) factors in h^'^\ Note 
that here we are treating the low energy Yukawa couplings as random variables, which are related to 
the couplings at higher energy scales via renormalization group equations. However, we expect that 
the effect of this renormalization group running can essentially be absorbed into a redefinition of the 
effective scale vb-l, hence leaving the results presented in the following unchanged. 
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chosen to be uniformly distributed in [0,27r). In the following, we shall refer to those 
sets of coefficients which are consistent with the experimental constraints in Eq. flTTl) 
as hits. 

In a preliminary run, we consider the neutrino mixing matrix U, with the effective 
scale vb^l = V~'^""^b-l/ sin^ /3 treated as random variable in the interval [lO""*^/^, 10^^^] x 
10^^ GeV. We find that the percentage of hits strongly peaks at vb-l — 1 x 10^^ GeV. 
This is interesting for two reasons. Firstly, it implies that given < a < 1, the high 
seesaw scale lies in the range 3 x 10^^ GeV < vb-l/ sin^ /3 < 1 x 10^^ GeV. Note that 
the upper part of this mass range is close the GUT scale, which is important for recent 
work on the connection of leptogenesis, gravitino dark matter and hybrid infiation [T^ . 
Secondly, this result allows us to fix the parameter vb-l in the following computations 
without introducing a significant bias. 

In the main run, for fixed vb-l, we include the mixing matrix Vl of the charged 
leptons to compute the full PMNS matrix. We require the mass ratios of the charged 
leptons to fulfill the experimental constraints up to an accuracy of 5% and allow for 
1 < tan /3 < 60 to achieve the correct normalization of the charged lepton mass spec- 
trum. Finally, imposing the 3a constraints on the two large mixing angles of the full 
PMNS matrix, we find parameter sets of 0{1) factors which yield mass matrices fulfill- 
ing the constraints in Eq. ( TTTj) . Our final results are based on roughly 20 000 such hits. 
For each hit we calculate the observables in the neutrino sector as well as parameters 
relevant for leptogenesis. The resulting distributions are discussed below. 

Statistical analysis 

In our theoretical setup the relative frequency with which we encounter a certain value 
for an observable might indicate the probability that this value is actually realized 
within the large class of fiavour models under study. In the following we shall therefore 
treat the distributions for the various observables as probability densities for continuous 
random variables. That is, our predictions for the respective observables represent best- 
guess estimates according to a probabilistic interpretation of the relative frequencies. 

For each observable we would like to deduce measures for its central tendency and 
statistical dispersion from the respective probability distribution. Unfortunately, it is 
infeasible to fit all obtained distributions with one common template distribution. Such 
a procedure would lack a clear statistical justification, and it also appears impractical 
as the distributions that we obtain differ substantially in their shapes. We therefore 
choose a different approach. We consider the median of a distribution as its centre 
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and we use the 68 % 'confidence' interval around it as a measure for its spread. Of 
course, this range of the confidence interval is reminiscent of the la range of a normal 
distribution. 

More precisely, for an observable x with probability density / we will summarize 
its central tendency and variability in the following form p^, 

X = f , A± = x± — X . (12) 

Here, a;_ and x+ denote the 16 %- and 84 %-quantiles with respect to the density 
function /. The central value x is the median of / and thus corresponds to its 50 %- 
quantile. All three values of x can be calculated from the quantile function Q, 

Q{p) = inf {x e [xmin, a^max] : P < F{x)} , F{x) = dt f{t) , (13) 

where F stands for the cumulative distribution function of x. We then have: 

x_ = g(o.i6), x = g(o.5o), x+ = g(o.84). (14) 

Intuitively, the intervals from Xmm to x_, x, and x+ respectively correspond to the x 
ranges into which 16 %, 50 % or 84 % of all hits fall. This is also illustrated in the 
histogram for sin^ 26'i3 in Fig. 1. Moreover, we have included vertical lines into each 
plot to indicate the respective positions of x_, x, and x+. 

In our case the median is a particularly useful measure of location. First of all, it is 
resistant against outliers and hence an appropriate statistic for such skewed distribu- 
tions as we observe them. But more importantly, the average absolute deviation from 
the median is minimal in comparison to any other reference point. The median is thus 
the best guess for the outcome of a measurement if one is interested in being as close as 
possible to the actual result, irrespective of the sign of the error. On the technical side 
the definition of the median fits nicely together with our method of assessing statistical 
dispersion. The 68 % confidence interval as introduced above is just constructed in such 
a way that equal numbers of hits lie in the intervals from a;_ to x and from x to x^, 
respectively. In this sense, our confidence interval represents a symmetric error with 
respect to the median. 

As a test of the robustness of our results, we checked the dependence of our distri- 
butions on the precise choice of the experimental error intervals. The results presented 
here proved insensitive to these variations. For definiteness, we therefore stick to the 
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Figure 1: Neutrino mixing angles 613 and 6*23. The vertical lines denote the position of the median (solid 
line) and the boundaries of the 68% confidence region (dashed lines) of the respective distribution. 



3a intervals. We also checked the effect of taking the random 0{1) factors to be dis- 
tributed uniformly on a linear instead of a logarithmic scale. Again, the results proved 
to be robust. 



4 Observables and results 

Mass hierarchy 

An important open question which could help unravel the flavour structure of the 
neutrino sector is the mass hierarchy. Since the sign of Aml^j^ is not yet known, we 
cannot differentiate with current experimental data between a normal hierarchy with 
one heavy and two light neutrino mass eigenstates and an inverted hierarchy, which has 
two heavy and one light neutrino mass eigenstate. Measuring the Mikheyev-Smirnov- 
Wolfenstein (MSW) effect of the earth could resolve this ambiguity. 

With the procedure described above, all hits match the structure of the normal 
hierarchy and there are no examples with inverted hierarchy. It is however notable 
that imposing the structure of the neutrino mass matrix given by Eq. ([7]) alone does 
not exclude the inverted mass hierarchy. Only additionally imposing the measured 
bounds on the mixing angles rejects this possibility. 

Mixing angles 

The mixing in the lepton sector is described by the matrix ?7pmns given in Eq. flTU]) . 
Of the three angles, two are only bounded from one side by experiment: for the largest 
mixing angle 623 there exists a lower bound, whereas the smallest mixing angle ^13 
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Figure 2: Lightest neutrino mass nii and effective neutrino mass in tritium decay m^. Vertical lines 
and shadings as in Fig. [TJ 



is SO far only bounded from above. Recent results from T2K [8J, Minos [9] and the 
preliminary result of Double Chooz [TU] point to a value of 6*13 just below the current 
experimental bound. The respective best fit points, assuming a normal hierarchy, are 
sin2 2ei3 = 0.11 (T2K), 2 sin^ ^23 sin^ 2^13 = 0.041 (MINOS) and sin^ 2^13 = 0.085 
(Double Chooz). The 90% and 68% confidence regions respectively read 

0.03 < sin^ 2^13 < 0.28 T2K, 90 % CL, Sep = 0, 

2 sin^ 023 sin^ 2^13 < 0.12 MINOS, 90 % CL, 6cp = 0, (15) 

0.01 < sin^ 2^13 < 0.16 Double Chooz, 68 % CL. 

With the procedure described above, we find sharp predictions for the smallest and 
the largest mixing angle within the current experimental bounds, 

sin^ 2^13 = Omt'om , sin^ 2^23 = 0.97^°;°^ ; (16) 

the corresponding distributions are shown in Fig.[TJ These results are quite remarkable: 
the atmospheric mixing angle points to maximal mixing, while the rather large value 
for 6*13 is consistent with the recent T2K, Minos and Double Chooz results. 

In our Monte-Carlo study we observe that the dominant contribution to the strong 
mixing in the lepton sector is primarily due to the neutrino mass matrix m^. The 
numerical results are not much affected by including the charged lepton mixing matrix 
Vl. The PMNS matrix is thus approximately given by the matrix U which diagonalizes 
the light neutrino mass matrix m^. 
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Figure 3: Effective mass in neutrinoless double-beta decay moup/s and Majorana phase a2i- Vertical 
lines and shadings as in Fig. [TJ 



Absolute mass scale 

The absolute neutrino mass scale is a crucial ingredient for the study of neutrinoless 
double-beta decay and leptogenesis. Although inaccessible in neutrino oscillation ex- 
periments, different experimental setups have succeeded in constraining this mass scale. 
Cosmological observations of the fluctuations in the cosmic microwave background, of 
the density fluctuations in the galaxy distribution and of the Lyman-a forest yield a 
constraint for the sum of the light neutrino masses, weighted by the number of spin 
degrees of freedom per Majorana neutrino, = 2, [TT] 

mtot = J]ym, <0.5eV. (17) 

V 

The Planck satellite is expected to be sensitive to values of mtot as low as roughly 
0.1 eV [13] • A further constraint arises from measuring the /3-spectrum in tritium 
decay experiments. The current bound [11] is 

= 5Z \^^PMNs)ei\^m\ < 4 eV^ (18) 

i 

By comparison, the KATRIN experiment, which will start taking data soon, aims at 
reaching a sensitivity of 0.04 eV^ [T^]. Finally, the neutrino mass scale can also be 
probed by neutrinoless double-beta decay. The relevant effective mass is 

m^vfifi = I y^jUpMNsfeimil ■ (19) 

i 

Here, Ref. [TB] claims a value of 0.11 — 0.56 eV. Dedicated experiments, such as 
GERDA [T7j with a design sensitivity of 0.09 — 0.20 eV, are on the way. Note that 
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Figure 4: Effective neutrino mass of the first generation nii and CP violation parameter ei. Vertical 
lines and shadings as in Fig. [TJ 



does not only depend on the absolute neutrino mass scale and the mixing angles, 
but also on the phases {a^i — 26) and 0:21 in the PMNS matrix. 

We find sharp predictions for the neutrino mass parameters discussed above. The 
lightest neutrino, z/i, is found to be quite light, cf. Fig. O 

mi = 2.21};^ X IQ-^eV, (20) 

hence favouring a relatively low neutrino mass scale beyond the reach of current and 
upcoming experiments. More precisely, we find for the neutrino mass parameters dis- 
cussed above: 

mtot = 6.01°:;^ X 10-2 eV, = 8.6tli x 10"^ eV, = 1.5+°| x 10~^ eV. (21) 

CP-violating phases 

The small value of the mass parameter measured in neutrinoless double-beta decay, 
^Quppi is due to the relative minus sign between the mi and m2 terms in Eq. (fT9|) . 
caused by a strong peak of the value for the Majorana phase 021 at vr, 

— = im- (22) 

vr 

This is depicted in Fig. El An analytic analysis of how this phenomena arises from the 
structure of the neutrino mass matrix, cf. Eq. ([7j), is presented in Appendix |X1 For the 
other Majorana phase asi and the Dirac phase 5 we find no such distinct behaviour 
but approximately flat distributions. 
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Leptogenesis parameters 

Finally, leptogenesis [IB] links the low energy neutrino physics to the high energy physics 
of the early universe. The parameters that capture this connection are the effective 
neutrino mass of the first generation rhi and the CP violation parameter ei [19] . 



j — 2,3 

with F{x) = y/x (in^^ + and Mj denoting the masses of the heavy neutrinos. 
Here, rhi determines the coupling strength of the lightest of the heavy neutrinos to 
the thermal bath and thus controls the significance of wash-out effects. It is bounded 
from below by the lightest neutrino mass nii. The absolute value of the CP violation 
parameter ei is bounded from above by [20] 



Sir v%y^ sm^ (3 Vsm /? / V 10 GeV / 

With the procedure described above, we find 

= 4.0+1^ X 10-2 eV , ^ = 0.25+°;28 , (25) 

and hence a clear preference for the strong wash-out regime Notice that there 

typically is a hierarchy between rhi and nii of about one order of magnitude. The 
relative frequency of the CP violation parameter ^i peaks close to the upper bound 
f^max, with the majority of the hits lying within one order of magnitude or less below 
^max5 cf- Fig- in This justifies the use of ^max when estimating the produced lepton 
asymmetry in leptogenesis. Here, in the discussion of ei, we assumed hierarchical 
heavy neutrinos, M2,3 ^ Mi. 

Theoretical versus experimental input 

The results of this section are obtained by combining two conceptually different inputs, 
on the one hand the hierarchy structure of the neutrino mass matrix given by Eq. ([T]) 
and on the other hand the experimentally measured constraints listed in Eq. flTT]) . In 
general, the distributions presented above really arise from the interplay between both 
of these ingredients. For example, the hierarchy structure alone does not favour a large 
solar mixing angle 6*12 and the ratio Am^Qj/Am^^j^ tends to be too large (cf. [^[22]). 
This discrepancy is eased by generating the random coefficients in Eq. ([T]) via the 
seesaw mechanism. Imposing the experimental constraints finally singles out the subset 
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of parameter sets used for the distributions presented above. As another example, 
consider the smallest mixing angle 6'i3 and the lightest neutrino mass eigenstate mi. In 
these cases, the hierarchy structure of the neutrino mass matrix automatically implies 
small values, similar to those shown in the distributions above. However, the exact 
distributions including the precise position of the peaks only arise after implementing 
the experimental constraints. A notable exception to this scheme is the Majorana phase 
0:21. Here the peak at 0:21 = vr is a result of the hierarchy structure of the neutrino 
matrix alone, as demonstrated in Appendix [Al 

5 Discussion and outlook 

In summary, we find that starting from a flavour symmetry which accounts for the 
measured quark and lepton mass hierarchies and large neutrino mixing, the present 
knowledge of neutrino parameters strongly constrains the yet unknown observables, in 
particular the smallest mixing angle ^13, the smallest neutrino mass mi, and the Majo- 
rana phase 0^21 . This statement is based on a Monte-Carlo study: Treating unspecified 
(9(1) parameters of the considered Froggatt-Nielsen model as random variables, the 
observables of interest are sharply peaked around certain central values. 

We expect that these results hold beyond Froggatt-Nielsen flavour models. An 
obvious example are extradimensional models which lead to the same type of light neu- 
trino mass matrix (cf. [23]). On the other hand, quark-lepton mass hierarchies and the 
presently known neutrino observables cannot determine the remaining observables in 
a model-independent way. This is illustrated by the fact that our present knowledge 
about quark and lepton masses and mixings is still consistent with ^13 ~ as well as 
with an inverted neutrino mass hierarchy (cf. [24j). As a consequence, further measure- 
ments of neutrino parameters will be able to falsify certain patterns of flavour mixing 
and thereby provide valuable guidance for the theoretical origin of quark and lepton 
mass matrices. 
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A Analytic derivation of the Majorana phase 0^21 



The complex phases of the 0{1) coefficients in the neutrino mass matrix rrii, and the 
lepton mass matrix rrie are randomly distributed. One would thus naively expect that 
also the Majorana phases 0:21 and a^i in the PMNS matrix can take arbitrary values. 
By contrast, the distribution of values for 021 that we obtain from our numerical Monte- 
Carlo study, cf. Fig. |3l clearly features a prominent peak at 0:21 = vr. In this appendix 
we shall demonstrate by means of a simplified example how the structure of the neutrino 
mass matrix m,^ may partly fix the phases of the corresponding mixing matrix U. 
Consider the following simplified Majorana mass matrix mj, for the light neutrinos. 





( 






nil, = V 
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(26) 



where (/? is an arbitrary complex phase between and 27r. For simplicity, let us neglect 
any effects on the mixing matrix U from the diagonalization of irie. That is, we define 
U such that U^nii^U = diag(mj), with denoting the eigenvalues of mlm^, 



rf sin^ (^/2) 2Tr7(5 + 3cos(y^))'/' + O [t]") 
^=4(l + r^^ [l-sinM^/2)])+0(r^^) . 



(27) 



Notice that the first two mass eigenvalues are nearly degenerate. This is a consequence 
of the particular hierarchy pattern of the matrix nii, which originally stems from the 
equal flavour charges of the and 63 multiplets. The relative sign of the O {rf) 
contributions to m{ and eventually shows up again in entries of U , for instance, 



U 



11,12 



2 (5 + 3 cos (^))^/^ 



--Arg [tz] 



(2J 



with z = 1 — cos {{p) — 2i sin (</?) . The phase 0:21 = 2 (Arg [f/12/f/ii] mod tt) in the matrix 
U represents the analog of the Majorana phase 021 in the PMNS matrix, cf. Eq. flTOj) . 
According to our explicit results for Uu and U12 it is independent of the arbitrary phase 
ip to leading order in r], 



^21 ~ 2 Arg 



exp 



1 X 

--Arg [+z\ + -Arg [-z] 



mod 



7r 



TT 



(29) 
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In a similar way we may determine the phase analogous to the Majorana phase a^i. 
However, due to the hierarchy between the mass eigenvalues mi and m^, the first and 
third column of the matrix U differ significantly from each other, thus leading to a 
phase that depends on ip at all orders of t]. 

Including corrections to all orders in rj and scanning over the phase (p numerically 
shows that the maximal possible deviation of 021 from vr is, in fact, of order 77^. Adding 
more complex phases to the matrix in Eq. (1261) gradually smears out the peak 
in the distribution of 021 values. The distribution that is reached in the case of six 
different phases is already very similar to the one in Fig. [31 We conclude that despite 
the need for corrections the rough picture sketched in this appendix remains valid: The 
hierarchy pattern of the neutrino mass matrix directly implies that 021 tends to be 
close to 021 = TT. 
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